Preparation of a Single Photon in Time-bin Entangled States via Photon Parametric 

Interaction 
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A novel method for preparation of a single photon in temporally-delocalized entangled modes is 
proposed and analyzed. We show that two single-photon pulses propagating in a driven nonabsorb- 
ing medium with different group velocities are temporally split under parametric interaction into 
well-separated pulses. In consequence, the single-photon "time-bin-entangled" states are generated 
with a programmable entanglement easily controlled by driving field intensity. The experimental 
study of nonclassical features and nonlocality of generated states by means of balanced homodyne 
tomography is discussed. 
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Entanglement and nonlocal correlations, besides their 
fundamental importance in the modern interpretation of 
quantum phenomena Q, are the basic concepts for 
realization of quantum information procedures The 
entanglement between matter and light states is an es- 
sential element of quantum repeaters |^ , the intermedi- 
ate memory nodes in quantum communication network 
aimed at preventing the photon attenuation over long 
distances. The two-photon entanglement is a crucial in- 
gredient for quantum cryptography 0, |(| , quantum tele- 
portation 0, 0, 0J , and entanglement swapping 0, ^| , 
which have been successfully realized during the last 
decade by utilizing two approaches, one based on contin- 
uous quadrature variables and the other using the polar- 
ization variables of quantized electromagnetic field [l2| . 
An essential step has been recently made in this direc- 
tion by implementing robust sources producing the pairs 
of photons which are entangled in well-separated tempo- 
ral modes (time-bins) ^ has been shown 0, 0] 
that this type of entanglement, in contrast to other ones, 
can be transferred over significant large distances without 
appreciable losses that makes it much preferable for long- 
distance applications. From the fundamental viewpoint, 
of special interest is a single photon delocalized into two 
distinct spatial ; 15j or temporal modes, for which case 
the nonlocality of quantum correlations is directly evi- 
dent from the violation of Bell's inequality formulated 
for the two-mode Wigner function ^(| . This was verified 
experimentally by performing the homodyne detection of 
delocalized single-photon Fock states and reconstructing 
the corresponding Wigner function from homodyne data 
[IT], Ell EH- To date two approaches have been developed 
for preparation of a single-photon in two distinct tempo- 
ral modes. In first one a time-bin qubit is created with 
the help of linear optics by passing a short pulse through 
Mach-Zehnder interferometer with different-length arms 
[l3|. The second approach is based on conditional mea- 



surement on quantum system of entangled signal-idler 
pairs generated via spontaneous parametric down con- 
version (SPDC) of successive pump pulses in a nonlin- 
ear crystal, when a detection of one idler photon tightly 
projects the signal field into a single-photon state coher- 
ently delocalized over two temporal modes [l9| . 

In this paper we demonstrate a novel method for dy- 
namical preparation of time-bin qubit. The basic idea 
is to create a parametric interaction between two single- 
photon pulses, which propagate in a driven medium with- 
out absorption and with slow, but different group veloci- 
ties. Then, due to the cyclic parametric conversion of the 
fields and the group delay, each pulse experiences a tem- 
poral splitting into well-separated subpulses. Moreover, 
since the process is completely coherent, at the output of 
the medium the time-delocalized single-photon states are 
formed. A remarkable feature of our scheme is the ability 
to produce two and more output temporally-entangled 
modes. Another important advantage is a generation in 
a simple manner any desired entanglement by controlling 
the driving field intensity. 

We consider an ensemble of A-type cold atoms with 
level configuration as in Fig.l. Two quantum fields 
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co-propagate along the z axis and interact with the atoms 
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FIG. 1: (a) Level scheme of atoms interacting with quantum 
fields Ei,3 and classical rf driving field of Rabi frequency Q. 
(b) Geometry of fields propagation. 
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on the transitions | ) — >| 1 ) and | ) — >| 2 ), 
respectively, while the electric-dipole forbidden transi- 
tion | 1 ) — >| 2 ) is driven by a classical and constant 
radio-frequency field (rf)with Rabi frequency f2 inducing 
a magnetic dipole or an electric quadruple transition be- 
tween the two upper levels. Here the electric fields are 
expressed in terms of the operators £i{z,t) obeying the 
commutation relations 



[£i{z,t),Eftz ,t)]=LS ij S(z-z) 



(1) 



where L is the length of the medium. We describe the 

N z 

latter using atomic operators & a p(z, t) = -A- J] | a )j(/3 | 

* i=i 

averaged over the volume containing many atoms N z — 
» 1 around position z, where N is the total number 
of the atoms. In the rotating wave picture the interaction 
Hamiltonian is given by 



Tr:Po2 = -r 2 /S 2 + igiiivm + zfTpoie lAkz - ig\£\p\% 
at 

(6) 
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P12 = + r 2 )/5i2 - igi£*po2 + ^2^2/Sio (7) 



Here 2Ti 2 are the decay rates of the excited states | 1 
) and | 2 ) and Ak = k 2 — fci — feii is the wave-vector 
mismatch. 

Further, we assume that the phase-matching condition 
Afc = is fulfilled in the medium. Then, the solution to 
Eqs.(5-7) to the first order in is readily found to be 



(8) 
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Here fcy = kd& z is the projection of the wave-vector of 

the driving field on the z axis, g a = PaoyJ^^ ^ s the 
atom-field coupling constants with p a p being the dipole 
matrix element on the transition | a ) — *| /3 ) , and V is 
the quantization volume taken to be equal to interaction 
volume. For simplicity, we discuss the case of exactly 
resonant interaction with all fields and, therefore, put in 
Eq.(l) the frequency detunings equal to zero, neglecting 
so the Doppler broadening, which in a cold atomic sam- 
ple is smaller than all relaxation rates. Then, using the 
slowly varying envelope approximation, the propagation 
equations for the quantum field operators take the form: 
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where Fi(z,t) are the commutator preserving Langevin 
operators, whose explicit form is not of interest here. 

In the weak-field (single-photon) limit, the equation of 
atomic coherences poi = <5oie~ 4 ' CiZ , i = 1,2 and p\i — 
di 2 e- li - k2 - k ^ z are treated pcrturbatively in £\ t2 - In first 
order only <t o — 1 is different from zero and for these 
equations we get: 
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Poi = -FiPoi+igi£iVoo + iQ* Po2e 



-ig 2 £2P2i (5) 



p 2 = /3oi(1^2), D = Q 2 +T 2 



(9) 



where, for simplicity, the optical decay rates are taken to 
be the same: Ti = T 2 = T. The first terms in right hand 
side (RHS) of Eqs.(8,9) are responsible for linear absorp- 
tion of quantum fields and define the field absorption co- 
efficients ki = 9i ™ upon substituting these expressions 
into Eqs.(3,4). Here the condition of electromagnetically 
induced transparency (EIT, refs.^ E3) ft > Tx. 2 is 
assumed to be satisfied for both transitions with weak- 
field coupling. The second terms in RHS of Eqs.(8,9) 
represent the dispersion contribution to the group veloc- 
ities of the pulses, while the two rest terms describe the 
parametric interaction between the fields. We require 
the photon absorption to be strongly reduced by impos- 
ing the condition kiL <C 1. Another limitation follows 
from AlueitT > 1 indicating that the initial spectrum 
of quantum fields is contained within the EIT window 
Aujeit — 7^ 775 where T is a duration of weak- 
field pulses, a — MuL is optical depth, a = ^:A 2 is 
resonant absorption cross-section and Af is the atomic 

number density. Finally, the length of the pulses has to 

o 2 

fit the length of the medium: Tvi < L with Vi = 

being the group velocity of the i-th field. Taking into 
account that kiL ~ T 2 a/£l 2 , this set of limitations yields 



r 2 



>> a 



and 
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(10) 



It is worth noting that upon satisfying the conditions 
(10), the dominant contribution to the parametric cou- 
pling between the photons is the third term in RHS 
of Eq.(8,9), because in this case the last term becomes 
strongly suppressed by the factor il 2 T/T >> 1. 

It is useful at this point to consider numerical estima- 
tions. The sample is chosen to be 85 Rb vapor with the 
ground state 5Si/ 2 (F g = 3) and exited states 5P 3 / 2 (F e = 
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2), bP 3 / 2 (F e — 3) of Rb atom as the states | ) and 
| 1 ),| 2 ) in Fig.l , respectively, using the following 
parameters: light wavelength A ~ 0.8fim, T = 2tt * 3 
MHz, atomic density Af ~ 10 12 cm~ 3 in a trap of length 
L ~ 100 /xm, SI ~ 10r, and the input pulse duration T ~ 
2-j-3ns. In this case a ~ 16, i>2 ~ 10 4 m/s, ui ~ 0.3v2, 
and fcji ^ 0.2. All of the parameters we use in our calcu- 
lations appear to be within experimental reach, includ- 
ing initial single-photon wave packets with a duration of 
several ns satisfying the narrow-line limitation discussed 
above. The standard method for producing single pho- 
tons based on SPDC in nonlinear crystals does not fit our 
purpose due to too broad linewidth (~ lOnm) of gener- 
ated light. Recently, a source of narrow-bandwidth, fre- 
quency tunable single photons with properties allowing 
to excite the narrow atomic resonances has been created 

mm. 

Then, taking into account that in the absence of pho- 
ton losses the noise operators Fi in Eqs.(5) give no con- 
tribution, the simple propagation equations for the field 
operators are finally obtained: 
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£x(z,t) = -i(3£ 2 
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(11) 



(12) 



where = gig 2 N/cQ is the parametric coupling con- 
stant. It is easy to check that these equations preserve 
the commutation relations (1). Note that for the pa- 
rameters above, the parametric interaction between the 
photons is highly strong /3L ~ 3. 

The formal solution of Eqs. (11,12) for field operators 
in the region ^ z L is written as 



£i{z,t) =£;(0, t—z/Vi)+ j dx{£ t (Q, t-z/vi 

o 
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-(z-x)) J (V)}, (13) 



where i,j = 1,2 and j ^ i. The Bessel function Jo(tp) 
depends on z via ip — 2/3y x(z — x) , Avij = Vi — Vj is 
the difference of group velocities. 

We are interested in dynamics of input state | ipi n ) = 
I li )<E> | O2 ) containing one photon in wi field. The 
similar results are clearly obtained in the case of one in- 
put photon at lo 2 frequency. We assume that initially the 
lo\ pulse is located around z = with a given temporal 
profile fx (i) : 

(0 I £x(0,t) I ipin) = (0 I £x(0,t) I lr) = fx(t) (14) 

The intensities of quantum fields at any time are given 
by 



Using Eqs. (13-15) and recalling that (0 | £2(0,*) | ipin) = 
0, we calculate numerically and show in Fig. 2 the 
output pulse shapes at z — L for the three values of 
£1 and for the case of Gaussian input (at z — 0) pulse 
fi(t) = cxp[— 2t 2 /T 2 } . For one-photon initial state, 
as is the case here, one can clearly see that the sec- 
ond field is not practically generated, thus demonstrat- 
ing that our scheme enables to prepare a single-photon 
in a pure temporally-delocalized state with an efficiency 
~ 100%. Moreover, depending on the driving field in- 
tensity, a different degree of initial pulse splitting and, 
hence, of entanglement is attainable. It is evident also 
that the total number of photons which is determined by 
the areas of the corresponding peaks is conserved upon 
propagation through the medium. Besides, in this case 
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FIG. 2: The numerical solution to Eqs. (15) at the output of 
the medium z — L for three values of fi. In these figures, 
solid curves represent the results for u)\ pulse, dashed curves 
show the UJ2 field generated in the medium, and dotted lines 
correspond to initial Gaussian pulse at uj\ frequency with T = 
2ns propagating in the medium in the absence of parametric 
interaction (3 = 0. For the rest of parameters see the text. 
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only two well-separated output temporal modes at oj\ 
frequency are produced, where due to v 2 > V\ a newly 
generated component is advanced compared to the signal 
pulse. This separation depends on the relative velocity of 
quantum fields, the larger the ratio V2/V1 the larger the 
group delay and the larger the output pulses separation. 
On the contrary, in the limit of equal group velocities the 
propagating pulses experience no splitting, as it can be 
easily seen from Eqs.(13), which in this case are reduced 
to 

£i(z, t) = £i(0, t)cos{/3z) - i£j(Q, r)sm(/3z) (16) 

where r = t — z/v\, j 7^ i. 

The system displays, however, a much richer dynamics 
in the case of input state | ipi„ ) = \ li )® | I2 ) consisting 
of one-photon wave packets at both frequencies oj\ and 
uj 2 . These results will be published elsewhere. Here we 
note only that in this case two multi-time-bin qubits at 
different frequencies io\ and lo 2 are generated, being at 
the same time strongly correlated with each other. This 
is evident also from the particular result of Eq.(16). 

The single-photon states are completely described by 
their Wigner function, whose remarkable property is that 
it takes negative values at the origin of phase space for the 
complex field amplitude. The negativity of the Wigner 
function is the ultimate signature of non-classical nature 
of these states. Besides, the nonlocality of quantum cor- 
relations between the two temporal modes directly fol- 
lows from the violation of Bell's inequality — 2 ^ B ^ 2 
predicted by local theories [0| • Here the combination B 
has the form: 

B = —[W(0, 0) + W(a u 0) + W(0, a 2 ) - W(a 1 ,a 2 )} 
where 

W( ai ,a 2 ) = \[2 \ai+a 2 | 2 _i] e -2M 2 -2| Q2 | 2 (17) 



is the Wigner function of two temporal modes calculated 
for the values of complex amplitudes ai = xi + iyi with 
Xi and y i} i = 1,2, being the quadratures of the i-th 
mode. In Eq.(16) we have supposed zero relative phase 
between superposition amplitudes of the two modes. In 
our case, for experimental verification of Bell's theorem, 
the Wigner function can be derived from the data of 
homodyne detection of quantum fields, when the sig- 
nals at the detectors are measured at two different times 
matched to the time separation between two toi output 
pulses obtained in Fig. 2. 



In conclusion, we have studied a highly efficient scheme 
for dynamic preparation of a single photon in distinct 
temporal modes, employing strong parametric interac- 
tion between two single-photon pulses, under the condi- 
tions of EIT, and their group delay. We have found the 
solution of propagation equations for the field operators 
depending on the propagation distance in terms of the 
Bessel function, the oscillatory character of which is just 
responsible for pulse temporal splitting. We have shown 
the ability of our scheme to achieve an arbitrary entan- 
glement by adjusting the driving field intensity, while the 
separation between the time bins can be controlled by us- 
ing the different atomic-level configurations to obtain the 
different ratio of group velocities of quantum fields. Sub- 
sequent papers will discuss the more complicated case of 
two input single-photon pulses and will present the re- 
sults of detail numerical simulations. 
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